Terahertz Bloch oscillator with a modulated bias 
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Electrons performing Bloch oscillations in an energy band of a dc-biased superlattice in the pres- 
ence of weak dissipation can potentially generate THz fields at room temperature. The realization of 
such Bloch oscillator is a long-standing problem due to the instability of a homogeneous electric field 
in conditions of negative differential conductivity. We establish the theoretical feasibility of stable 
THz gain in a long superlattice device in which the bias is quasistatically modulated by microwave 
fields. The modulation waveforms must have at least two harmonics in their spectra. 

PACS numbers: 03.65.Sq, 73.21. Cd, 07.57.Hm, 72.30.+q 



Electrically biased semiconductor superlattice (SL), 
which operates in conditions of a single miniband trans- 
port regime, exhibits a static negative differential con- 
ductivity (NDC) due to Bloch oscillations Assuming 
a homogeneous electric field inside the nanostructure, the 
semiclassical theory predicts that SL in the NDC state 
can provide small-signal gain for all frequencies below 
the Bloch frequency [2j. If miniband electrons can per- 
form relatively many cycles of Bloch oscillations between 
scattering events, the Bloch gain profile is shaped as a fa- 
miliar dispersion curve (Fig. [I]:, blue curve). In general, 
the dispersive gain profile is a quantum dissipative phe- 
nomena which is not limited to the Bloch oscillating elec- 
trons. Similar profiles with a resonant crossover from loss 
to gain have been observed in optical transitions of quan- 
tum cascade lasers Q and microwave response of Joseph- 
son junctions Importantly, estimates for typical SLs 
predict a significant THz gain even at room temperature 
for frequencies in the vicinity of Bloch gain maximum. 
Thus, the principal scheme of CW Bloch oscillator would 
consist of a biased SL which is placed in a high-Q cavity 
tuned to a desirable THz frequency. However, even for 
moderate carrier densities a state with NDC is unstable 
against a development of space-charge instability which 
eventually results in the formation of high-field domains 
inside a long sample The formation of electric do- 
mains, which breaks the assumption of a homogeneous 
electric field used in Q , ultimately represents one of the 
major obstacles to the realization of Bloch oscillator. 

In short SLs the formation of electric domains can be 
suppressed. However, since a short sample can emit only 
a weak electromagnetic power a stack of such SLs is re- 
quired to observe the Bloch gain Q . In our research we 
focus on the feasibility of Bloch gain in a long SL. To 
circumvent the electric instability in a long SL it is nec- 
essary to modify the unstable Bloch gain profile in such 
a way that the sample exhibits a positive dynamical con- 
ductivity at low frequencies, including dc, whereas the 
high-frequency conductivity is still negative giving rise 
to the desired amplification @, S] • To reach this goal it 
was suggested to use SLs with narrow gaps, where addi- 



tionally to the intraband Bloch oscillations the interband 
Zener tunneling becomes important Q . In a recent pub- 
lication we theoretically showed that the desirable 
modification of Bloch gain profile still can be achieved 
within the standard single miniband transport regime. 
We utilized additional peaks in the voltage-current (VI) 
characteristics of SL, which arise under the action of a 
strong ac field if the Bloch frequency is close to the ac 
field frequency or its harmonics However, the need 
for the use of a strong THz field as ac pump limits the 
applied aspect of this suggestion. 

In this Letter, we show that the unstable Bloch gain 
profile can be made stable by means of a low-frequency 
modulation of dc bias applied to a long SL. At least two 
harmonics with a proper phase difference are necessary 
to be present in the modulation spectrum in order to 
suppress the destructive low-frequency instability. The 
degree to which the large magnitude of THz gain near 
the Bloch resonance may be preserved here depends on 
the modulation waveform. In this respect we find the 
pulse modulation optimal. Since the operation of Bloch 
oscillator requires that during at least some portion of 
the modulation period SL is biased to the NDC state, 
an accumulation of charge can arise. The requirement 
that the charge accumulation is weak imposes a lower 
limit on the modulation frequency, which, however, still 
belongs to the microwave frequency domain. Finally, our 
numerical simulations show that the effect of stable THz 
gain persists if boundary conditions in the SL device are 
taken into account. 

We consider a SL under the action of an electric field 
E(t) — E mo d{t) + Egcosujt, where E$cosujt is a weak 
probe field with the frequency u) fixed by an external 
circuit (external cavity) and E moc i(t) is a periodic mod- 
ulation of the bias (Fig. [T] (d)). We suppose that the 
modulation frequencies are incommensurate with ui. To 
understand the main idea of our suggestion let us first 
consider E mo d(t) in the forms of pulses as shown in Fig. [I] 

(a) and (b) . For both the asymmetric (a) and symmetric 

(b) waveforms the amplitude of modulation Eq is such 
that during each period T SL spends most of the time 
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FIG. 1: (a), (b) Two different waveforms of bias modulation 
E mo d{i) which result in the same dynamical conductivity of 
SL [Eq. {T}]. (c) Gain profiles of SL calculated using Eq. ([T]) 
for E = 4E cr and T /T =: 0, 0.05, 0.1, 0.2. (d) Sketch of the 
Bloch oscillator with a modulated bias. 



T — To in the NDC region of the VI characteristic and 
only for a short fraction of the period Tp/T <C 1 the 
bias is switched off. To calculate the gain profile one 
needs to know the real part of the dynamical conductiv- 
ity a r (u>) = Re[cr(w)]. Gain corresponds to a r (io) < 0. 
We assume that the switching is faster than both T and 
T—Tq and transients after the switching can be neglected 
T — To, To t (r is the intraband relaxation time) 
We find for the both sequences of pulses 
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are the dynamical conductivity of a biased SL and free- 
electron absorption in an unbiased SL, respectively. Here 
jdc(eE dc d) is the Esaki-Tsu characteristic [l|, 



jdc(eEdcd) = j p 
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is the peak current corresponding to the critical field 
E cr = h/edr, Ik{x) (k — 0,1) are the modified Bessel 
functions, iV is the density of electrons in the first mini- 
band of the width A, d is the SL period, T e is the temper- 
ature and ap = 2j p /E cr is the Drude conductivity of SL. 



Note that the finite difference form 14[ given by Eq. ([2]) 
is equivalent [l3} to the original Ktitorov et al. (KSS) 
formula Q for a single relaxation time. 

In the limit T /T — * 0, Eq. |T]) describes the usual 
Bloch gain profile (thin blue line in Fig. Q] (c)). With 
an increase of Tp/T SL spends more time in an unbi- 
ased state, therefore the gain at low frequencies decreases 
and eventually the dynamical conductivity becomes pos- 
itive, whereas the high-frequency gain is still preserved 
as shown in Fig. 1 (c). The physical origin of such be- 
havior is quite intuitive: At high frequencies around lj ~ 
uob — t^ 1 (ujb = eEpd/h >• r _1 is the Bloch frequency), 
the Bloch gain a^ ss (uj) w —ap/AujT dominates the free- 
electron absorption af ree (uj) ps ap/uj 2 T 2 . In contrast, 
at low frequencies lot <C 1, the free-electron absorption 
a free w ffQ overcomes the Bloch gain a^ ss —ap/uo B T 2 . 

The dynamic conductivity expressions Eqs. dTJ-© di- 
rectly follow from the general formula describing the dy- 
namic conductivity of SL in the case of arbitrary slow 
(quasistatic) modulation 
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where j d n ° d is the dc component of the current in SL 
1 f T 
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that adiabatically follows to the time-dependent varia- 
tions E mo d(t) via the Esaki-Tsu characteristic 0$. These 
formulas can be derived using two different techniques. In 
the first method, we employed the standard semiclassical 
approach based on an exact solution of Boltzmann trans- 
port equation for a single tight-binding miniband and a 
constant relaxation time [13(. Next, we also found that 
the same dependence of <j(lo) on the shape of VI char- 
acteristic can be obtained in Tucker theory of quantum 
tunneling in the presence of ac electric fields [14| . Both 
approaches result in Eq. ([5]) with a cumbersome expres- 
sion for the jX° d ' which can be further reduced [15( to 
the simple form of Eq. © in the limit of low-frequency 
modulations. 

For the pulse modulation the spectrum of E mo d(t) 
consists of many harmonics of the basic frequency Q(= 
2ir/T) and, as we have seen, the gain profile has the 
desirable shape with cr r (u> — > 0) > 0. Importantly, 
in the case of monochromatic quasistatic field no high- 
frequency gain exists if a r {uj — > 0) > [lj|. Therefore, 
we conclude that the monochromatic low-frequency mod- 
ulation cannot provide the required stabilization of the 
space-charge instability in the Bloch oscillator. Now we 
are going to show that a bichromatic modulation can al- 
ready mimic the kind of fast bias switching necessary for 
the stabilization. We consider a modulation 



E mo d(t) = Edc + Ei cos fit + E 2 cos(nQt 
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FIG. 2: (a) Esaki-Tsu characteristic (thin blue line) [Eq. @] 
and the characteristic of SL modified by the bichromatic field 
[Eq. ©] (thick black line) for Ei = 2.lE cr , E 2 = 1.7E cr , 
Qt — 0.05, n = 3, (j> = 0. (b) Dynamical conductivity 
[Eq. ([5])] as a function of u> (thick black) for the same bichro- 
matic modulation and Edc = 4E cr . The thin blue curve shows 
the usual Bloch gain profile [Eq. @] . 



where n is an integer and Qt <C 1. Fig. [5] demonstrates 
the VI characteristic and gain profile for the modulation 
with moderate field strengths £1,2- The amplitude E2 
of the harmonic component and the relative phase <fi are 
chosen to square up the field (|7|). For typical scattering 
time r = 200 fs the modulation frequency used in the 
figure is £l/2ir ss 40 GHz. This E mo d(t) creates, follow- 
ing Eq. ©, an additional peak in the VI characteristic 
of SL (Fig. Ufa)), which resembles the resonant struc- 
tures in VI characteristic induced by a monochromatic 
THz field This additional peak can be utilized in a 
similar way as in the earlier theoretical suggestion 10] to 
achieve THz gain at a stable operation point. Namely, we 
put the working point at the part of the peak with a posi- 
tive slope and then apply the difference formula (0 • The 
result is the dispersive gain profile with a positive dynam- 
ical conductivity at low frequencies as shown in Fig.[2jb). 
In contrast to the earlier work [10|, for the suppression 
of the electric instability in SL we use a bichromatic mi- 
crowave field [jt| instead of a monochromatic THz field. 
Importantly, the shape of the VI characteristic and mag- 
nitude of gain can be additionally controlled by a varia- 
tion of the relative phase <p. The magnitude of THz gain 
at bichromatic modulation is a bit lower in comparison 
with the case of pulse modulation (cf. Figs. and 

mc)). 

Apart from a r (u> — > 0) > 0, the conditions for space- 
charge control in microwave driven SL with NDC include 
also the requirement SIt^ > 1 [13], where Td is a char- 
acteristic time of domain formation. It means that an 
accumulation of charge during the period of modulation 
T is always weak 17, Hq . If the conditions flTd > 1 
and <r r (u> — >• 0) > are satisfied, VI characteristic of SL 
is modified by strong microwave fields according to the 
quasistatic theory prediction with an assumption of spa- 
tially homogeneous fields, i.e. following Eq. |6]). Exper- 
imental observation of such modifications in the case of 
monochromatic microwave field has been reported in [l9| . 
The time Td is of order of the dielectric relaxation time 



seo/\(Jdc(Edc)\ [20|, where (Jdc(Edc) is the dc differential 
conductivity, e is the relative permittivity of SL material 
and eo is the permittivity of vacuum. For the Esaki- 
Tsu characteristic (Eq. Q) mm{<Jdc{Edc)] — — co/8, and 
therefore we can present the condition of limited charge 
accumulation in the form 



r 2 < 80r, 
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(2j p ed/heoc) is the miniband plasma fre- 
2l|. Since ujh oc NA, for given electron con- 



centration N and miniband width A, the inequality ([8]) 
imposes a lower limit on the modulation frequency f2. Al- 
ternatively, for a fixed value of the modulation frequency 
fi, Eq. ([5]) requires the use of SLs for which the product 
NA is less then some critical value. 

The absorption (gain) a in units cm -1 is related to the 
scaled dynamical conductivity as 
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where c is the speed of light in vacuum. Comparing 
Eqs. (JHJ and © we see that the condition ([5]) imposes a 
restriction on the maximal value of a. For the parame- 
ters of bichromatic modulation considered in Fig. [2l we 
estimate a ~ —3.6 cm -1 at the gain resonance lot « 2.5. 
The corresponding value w^t 2 = 0.4 [Eq. ((8])] can be 
realized, for example, at room temperature in SL with 
d = 6 nm, A = 60 meV and N = 6 x 10 15 cm" 3 . The 
magnitude of gain can be increased by using the rectan- 
gular wave modulations. In this case, for the same basic 
modulation frequency fir = 0.05 we get a w —12.8 cm' 1 
at THz gain resonance shown in Fig. [IJc). It requires 
the use of ~ 10 ps dc pulses. However, a proof-of-the- 
principle demonstration, perhaps in the pump-probe con- 
figuration, can be performed with commercially available 
electric pulses of the duration ~ 0.1 ns. The restriction 
([8]) can probably be relaxed in the case of lateral surface 
SLs formed in two-dimensional electron systems. Be- 
cause the domain growth rate can be suppressed in 2D 
systems 23] , we can expect reasonably large magnitudes 
of THz gain even for a rather slow modulation of the bias. 

Degree of homogeneity of the electric field inside the 
SL also depends on the boundary conditions defined by 
the attached electric contacts. Since the THz probe 
field is weak and bias modulation is quasistatic, we can 
use the standard drift-diffusion model ItJ with typi- 
cal boundary conditions [24| . SL with a nonlinear con- 
ductivity, which follows from the Esaki-Tsu dependence 
(U]), is placed between two Ohmic buffer regions of the 
equal conductance ct& 24(. The ac voltage across the 



device is periodically varied as shown in Fig. QJa,b). 
We numerically calculated the degree of homogeneity 
as r = ([E max (t) - E min (t)} / E Q ) U where E max (t) and 
E m in{t) are the maximum and minimum of the electric 
field in the whole SL and the time- averaging (...)* is 
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FIG. 3: Regions of a spatially homogeneous electric field 
(dark, red online) in ab-Upi plane in the case of (a) asym- 
metric and (b) symmetric pulse voltages (see Fig. [T] (a) and 
(b), respectively) with T /T = 0.1. SL has the length 500d. 



performed over a long time interval of sustained field os- 
cillations. Fig. [3] shows the values of r in the o^-io^ plane 
for a very long SL. In the indicated dark regions of small 
r values the electric field inside the SL can be consid- 
ered as spatially homogeneous. We also verified that for 
these parameters spatially-averaged gain profiles [25| are 
practically indistinguishable from those calculated ear- 
lier (Fig. [He)) within the homogeneous field approxima- 
tion. If the buffer conductance coincides with the SL 
conductance (ab/&o » 0.06 in Fig. [3|), the upper bound- 
ary of the homogeneous field region in w^r 2 is even a bit 
larger than expected from Eq. ([8]). For mismatched con- 
ductances smaller values of to^ are required to achieve a 
desirable homogeneity. Interestingly, in the case of sym- 
metric modulation SL is less sensitive to the boundary 
effects (cf. (a) and (b) in Fig. [3]). We attribute it to 
the fact that the time-averaged current induced by the 
symmetric voltage is zero both in the buffers and in SL 
independently on the value of <7b and thus, on average, 
there should be no accumulation or depletion of charges 
at the boundaries. For shorter SLs we found that the 
difference between the effects of symmetric and asym- 
metric waveforms becomes much less prominent and, as 
expected, the field stays homogeneous for a much wider 
range of Ob- We underline that an observation of the 
Bloch gain certainly requires a gradual change of con- 
ductance between the active SL region and the metallic 
contacts. 

In summary, we have theoretically shown that a proper 
quasistatic modulation of the bias in Bloch oscillator re- 
sults in THz Bloch gain in conditions of positive differ- 
ential conductivity. Waveform of the modulation should 
include two distinct parts: Time interval of almost con- 
stant bias responsible for the THz gain due to the exci- 
tation of Bloch oscillations and a virtually unbiased in- 
terval during which the dominant free-carrier absorption 
effectively suppresses the undesirable low-frequency in- 
stability. Spectrum of the modulation must have at least 
two harmonics with a controlled difference of phases. 

Along with the semiconductor superlattices our results 
can be directly applied to such semi- and superconduct- 
ing microstructures as the dilute nitride alloys [26j and 



small Josephson junctions [27| . It has already been found 
that, despite of the very different physical origin, the 
small Josephson junction demonstrate in fact the Esaki- 
Tsu VI characteristics and ac response similar to the sin- 
gle miniband superlattice 271 ]. 
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